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Abstract 

We study a general field theory of a scalar field coupled to gravity througli a quadratic Gauss- 
Bonnet term ^{(j)) Rqb- '^^^ coupling function has the form (^{(p) = where n is a positive integer. 
In the absence of the Gauss-Bonnet term, the cosmological solutions for an empty universe and a 
universe dominated by the energy-momentum tensor of a scalar field are always characterized by 
the occurrence of a true cosmological singularity. By employing analytical and numerical methods, 
we show that, in the presence of the quadratic Gauss-Bonnet term, for the dual case of even n, 
the set of solutions of the classical equations of motion in a curved FRW background includes 
singularity-free cosmological solutions. The singular solutions are shown to be confined in a part 
of the phase space of the theory allowing the non-singular solutions to fill the rest of the space. We 
conjecture that the same theory with a general coupling function that satisfies certain criteria may 
lead to non-singular cosmological solutions. 
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INTRODUCTION 



Despite the successes of Einstein's theory of Gravitation at large distances, a Quantum 
Theory of Gravity, vahd at supersmall distances, requires a more general framework. At 
present, Superstring Theory |l|] seems to provide the most appealing framework for such a 
theory. Superstring Theory leads to the unification of gravity with the other fundamental 
forces. It also leads to important modifications of the standard cosmology, based on the 
Einstein action, at short distances of the order of the Planck length. Although the theory 
has not been fully developed to the point that a detailed cosmology could be constructed, 
a number of general conclusions can be drawn regarding new possibilities that distinguish 
string cosmology from the standard model. Modifications of gravity of stringy origin can 
be studied through the Superstring Effective Action corrected by incorporating loop and a' 
effects. The latter are associated with the contribution of the infinite tower of massive string 
modes, while the former are due to quantum loop effects. Although the full String Theory 
is approximated only in a perturbative sense by the Effective Action and this is expected 
to describe Physics only up to energies where quantum gravitational effects start becoming 
dominant, it is hoped that the loop-corrected Action captures many of the true features of the 
exact theory. The study of the loop-corrected Superstring Action has uncovered interesting 
possibilities not realized by the Einstein-Hilbert action such as the existence of novel 
stable dilatonic black holes that circumvent the "no hair" theorem in its restricted sense. 
There are of course alternative approaches to string cosmology. The absence of cosmological 
singularities in the presence of higher-curvature terms, in various numbers of spacetime 
dimensions, has been pointed out in the literature several times @] . Another approach is the 
pre-Big Bang scenario P which attempts to incorporate features of the exact theory such 
as duality symmetries. 

A remarkable property of the loop-corrected Superstring Effective Action in the presence 
of the dilaton and moduli fields is the existence of singularity-free solutions with flat initial 
asymptotics [0. These are linked to the gravitational terms with field- dependent coeffi- 
cients that are present. These solutions which avoid the initial singularity are possible for 
a definite sign of the corresponding trace anomaly for which the strong energy conditions 
related to the modulus energy-momentum tensor can be violated. They start from a fiat 
space-time in the infinite past, they pass through an inflationary period and they end up as a 
slowly expanding universe. A general field theory of a scalar field coupled to gravity through 
a quadratic Gauss-Bonnet term ^{(f))RQ^ has also been shown to possess singularity-free 
solutions in a spatially fiat FRW background under very mild assumptions on the coupling 
function ^(0) 0. In a subsequent paper by Easther and Maeda [|], the case of a closed 
FRW universe based on the loop-corrected Superstring Action was also shown to lead, via 
numerical methods, to such singularity-free solutions. 

In the article at hand, being inspired by the Superstring Effective Theory, we consider 
a generic theory with a scalar field coupled to gravitation through the higher-curvature 
quadratic Gauss-Bonnet term. For simplicity, we keep only the one-loop gravitational quan- 
tities that appear in the action functional of the Superstring Effective Theory. For the case 
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of a spherically symmetric background that we are going to consider, the other one-loop 
gravitational quantity RR = rj'^^'^'^ R"'^^^R^\^r vanishes identically. This model captures the 
essential features and was shown to possess non-singular solutions in the fiat case for a ge- 
neral class of coupling functions [0. Here, we extend the analysis of Ref. ||^ to the case of 
a curved universe, both open or closed. We develop a purely analytical argument and we 
manage to show that the singular cosmological solutions, characterized by a true singularity 
at finite time, are indeed present in the theory but they are confined in a certain part of the 
phase space of the theory. In this way, the non-singular cosmological solutions are summoned 
to fill the rest of the space. These results are radically different from those that follow from 
the same theory when the quadratic Gauss-Bonnet term is absent. In that case, as we will 
show, the singular cosmological solutions cover the whole phase space of the theory leaving 
no room for the existence of non-singular solutions. 

The structure of this article is as follows: In section 1, we derive the equations of motion 
for the scalar and gravitational fields in a curved FRW background. In section 2, we study 
the cosmological solutions of the theory when the Gauss-Bonnet term is absent. We consider 
both the cases of an empty universe and a universe dominated by the energy-momentum 
tensor of a scalar field. In section 3, we develop our analytical argument for the existence of 
non-singular cosmological solutions in the presence of the Gauss-Bonnet term. In section 4, 
a numerical analysis for a specific choice of the coupling function serves as an illuminating 
example for our theory. The last, short section is devoted to our conclusions. 



1 Equations of motion of the theory 

We consider the quadratic coupling of a scalar field with gravity through the Gauss- 
Bonnet term which is described by the action 

S = j d'x^g[^ + U,<pdy-^5mnlj}^ . (1.1) 

Note that ^((/)) is, for the time being, a general coupling function. The Gauss-Bonnet term 
is defined as 

T^Ib = R^.up.R'"'"' - ^Rp.R^" + R' (1.2) 

and S, which in Superstring Effective Theory is proportional to the trace anomaly of the 
theory, plays the role of a coupling parameter. 

The spacetime background assumes the standard spherical symmetric FRW form 

ds^ = dt" - e^-W I + {de^ + sin'e d^') | , (1.3) 

[ 1 — kr'' J 

where k = 0,±1 corresponding to the fiat, closed and open universe, respectively. 
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Making use of the above metric components and assuming further that the scalar field 
depends solely on the time coordinate t, the equations of motion take the form 



df 



+ 3(J)uj -24^(uj + uj'^) iuj'^ + ke 



-2uj\ 







;i-4) 



2 (d) + (1 + 8/cI;) + {Cj' + ke-'"^) (1 + 8/) + - 



3 {6j' + A;e-'") (1 + S/cI;) - ^ = , 

where = — (5^(0)/16. If we set x = 0, 2; = ci; and y = e~^'^, we obtain 

3 

x + 3zx + - 5^'{z + z^) (z^ + ky) = 

(i + z^) (4 - 26Cxz) + {z^ + ky) (2 - 5i"x^ - d^'x) + = 
{z^ + ky) (6 - 35^' xz) - = 



;i.5) 
;i.6) 



;i-7) 

1.8) 

;i.9) 



Rearranging eqs. ( |1.7| ) and ( |1.8| ), we obtain a new equation which contains only the time 
derivative of z 



dz 
z = — X 



2 - 5Cx^ + 36Czx) {z^ + ky) + x^ 
4 - 26Czx + I (5^)2(^2 + ky)^ 



[1.10) 



while, from the definition of y we are led to the following differential equation 

?/ = ^ a; = -2yz . 



[1.11) 



On the other hand, we may solve eq. (|1.9| ) as an algebraic equation and write x, the time 
derivative of the scalar field 0, as a function of z, y and 6^' in the following way 



X = --6^'z{z^ + ky) + s] 



5i'z (^2 + ky) 



n2 



+ 6 (^2 + ky) 



S = ±l 



(1.12) 



Note that the set of equations ( p..lO| )- (|l.ll| ) is characterized by an invariance under the 
simultaneous change of the signs of z and s. In order to clarify this point, we suppose that 
we have found a solution, for the choice s = +1, described by the set of equations 



dz 

'=d4>''^ 
dy 



(2 - 6Cxl + 35^'zx+) (^2 + ky) + x^^ 
4 - 26CZX+ + I (5^)2(^2 + ky)^ 



-2yz . 



[1.13) 
[l.U) 
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where x+ stands for the value of x that corresponds to the choice s = +1. Under the 
transformation z — —z, the set of equations ( p..l3| )-( p!.14| ) is replaced by a new one with 
a;_ in the place of X4., where x_ corresponds to the choice s = —1. This means that if a;+ 
corresponds to a solution (singular or not) of the equations of motion, then x_ corresponds 
to the same solution with the sign of z reversed. For this reason, we may keep fixed the sign 
of z, e.g. z > 0, during the analytical treatment of the problem. 



2 The (5 = case 



We first consider the case with 5 = 0, that is without the Gauss-Bonnet term. If we 
further assume that the scalar field takes on a constant value and set x = = in equations 
(|r|)-(|r|), we obtain 

+ ke-^"^ = Q , (2.1) 
[;j + u^ = {}. (2.2) 

For A; = 0, we obtain uj = const, which corresponds to a static universe with arbitrary 
radius. On the other hand, for k = +1, we must have u = const, and e~'^ = at the same 
time, which corresponds to a static universe with infinite radius. The only interesting case 
is the last one, k = —1, where we find that e'^^*-' = a{t) ~ t. This result corresponds to a 
linearly, eternally expanding universe with an initial singularity, at t = 0. Note that the rate 
of expansion is much larger than in the case of the "radiation" (a ~ t^^^) or the "matter" 
(a ~ t"^^^) epoch of the Standard Cosmological Model. This is due to the absence of any 
matter content of the universe capable of slowing down the expansion of the universe. 

Next, we allow the scalar field to evolve with time (0 7^ 0) while keeping the parameter 
S equal to zero. In this case, the set of equations (^TTO|)-(|LTTj) take the form 

dz 2 
z = — X = —3z — 2ky , 
d(p 

y = ^x = -2yz, (2.3) 

where 



X = s^Qiz"^ + ky) , s = ±1 . (2.4) 
We are going to study separately the cases of fiat and curved space : 

A) Flat Space (A; = 0) : The solution of (pl3|)-(p^) with respect to time t takes the 
form 

z = -Sz"^ z{t) = {c + 3ty^ , 
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y=-2yz y{t) = {c + 3ty^/^ , (2.5) 



4> = sV6z (j){t) = S]^-ln{c + 3t) + c' 

where c and c' are arbitrary constants. The resuh for the scale factor of the universe is 

e'^ = a(t) = (c + 3t)'/^ (2.6) 

which corresponds to an expanding universe with a true cosmological singularity at finite 
time. Note that, here, the rate of expansion of the universe is smaller than the corresponding 
ones during the two epochs of the Standard Cosmological Model. The sole reason for this 
result is the presence of the energy momentum tensor of the free scalar field on the right- 
hand side of the Einstein's equations which leads to the slowing down of the expansion of 
the universe in a more effective way than the energy momentum tensor of a perfect fluid. 



B) Curved Space (k = ±1) : In this case, the system (|2.3|) may be reduced to a single 
equation : 

dz dz , ^ 9 , , , 2A;-u , 

z = — y= — (-2yz) = -3z^ - 2ky =^ 2ydz - 3zdy = ^ dy . 2.7 
dy dy z 

If we multiply both sides by z/y^, we obtain 

^(^) = (-^) ^(^) ^z' + ky = ciy' , (2.8) 

where Ci is a positive constant. Substituting the above in the differential equation of y, we 
get the result 

dy dy 



y=—x = —sJ 6ci?/3 = -2y J y (ci?/2 - k) ^ (2.9) 



(2.10) 



Cik + C2 exp ■ 




2 ci C2 exp 1 




\ 



From the above expression as well as from eq. (|2.9|) , it is evident that there is a further 
invariance of the solutions under the interchange of the signs of s and 0. As a result, we may 
keep fixed the sign of s, e.g. s = +1, while allowing to take on both positive and negative 
values. 

A cosmological singularity is encountered when a{t) ^ or equivalently when y 00. 
From the expression ( |2.10|) , we conclude that, when k = +1, y goes to infinity for —>■ ±00 
while, for k = —1, a singular behaviour arises only for — > —00. Near the singularities, we 
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may evaluate the approximate expression of y which can be written in the following way 

For ^ -oo ^ ?/ ^ ;^e~v^* (for A; = ±1) , (2.11) 

For ^ +CX) y ~ — e^"^ (for A; = +1) . (2.12) 

2C2 

The corresponding expressions for z can be easily derived from eq. (|2.8|) and exist only when 
the solution for y exists as well. By making use of the differential equation for 0, eq. ( |2.4| ), 
we may deduce the dependence of the scalar field on time t near the singularities and, 
consequently, the expression of the scale factor of the universe in the same region. Then, we 
obtain 

a{t) ~ (c + 3v^t)^/^ (for A; = ±1) , (2.13) 

a{t) ~ (c' - ^^itf'^ (for k = +l) . (2.14) 

The above expressions describe also a universe with a true cosmological singularity at finite 
time. We note that for k = +1, that is for the case of a closed universe, there are always 
two branches of singular solutions with vanishing a{t). On the other hand, for the options 
k = 0,-1 which correspond to the cases of a flat and open universe, respectively, there is 
only one branch of singular solutions. This result is in perfect agreement with the singularity 
content of the Standard Cosmological Model. The open and fiat universe are characterized 
by only one singularity, the initial one, while in the case of a closed universe we encounter 
two cosmological singularities, the initial and the final one. 

It is also worth noting that, for the choice 6 = 0, the group of singular solutions found 
above covers the whole phase space of the theory leaving no space for the existence of non- 
singular solutions. The final singularity of the closed universe ( p.l4|) can be avoided only if 
we choose k = 0,-1. On the other hand, the initial singularity (|2.13| ) disappears only if we 
set C2 = 0. Then, we end up with the totally unrealistic case of a static universe with infinite 
radius. As a result, we conclude that, in the absence of the Gauss-Bonnet term, the only 
realistic cosmological solutions, that we may obtain in the framework of the theory ( |1 . 1| ) , 
contain, at least, one true singularity. 



3 The S^O case 

In this section, we are going to search for non-singular cosmological solutions in the presence 
of the quadratic Gauss-Bonnet term in the action functional of the theory. It will be useful 
for our analysis to search for violations of the energy conditions |^ that indicate absence of 
singularities. Assuming a perfect fluid form for the energy-momentum tensor of the system. 
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the energy and pressure are defined as : Too = Tu = —pgu- Using the equations of motion 
(|1.4|) - (|1.6|) , the energy conditions take the form 



, 2x B + 2Ax'^z^ -Q5i"x^z^{z'^ + ky) , , 

p + p= -2(c^-A;e-2-) = 2(A:?/ + 22) ^ v y; ^ ^3^^^ 

B 

p + 3p = _6 {Co" + 0)) = (8 - b^'x^) {ky + z^f , (3.2) 

B 

where 

B = x^ky + 5x^z^ -Ux^k^y^ -24:X^kyz^ -Ux^z^ + 

108 kyz'^ + 36k^y^ + 108 k'^ y^ z^ + 36 z^ (3.3) 

and where we have used eq. (1.5) in order to eliminate . In this form, we may easily 
prove that for k = 0, +1 the term B, being a polynomial with respect to with no real 
roots, is always positive definite. Thus the energy conditions can only be violated, leading 
to non-singular cosmological solutions, for S > 0. For k = —1, the analysis is much more 
complicated but it can be shown that the energy conditions are violated for both signs of 
6. However, by making use of numerical as well as analytical arguments, we may show that 
non-singular solutions arise only for 6 > 0, too. As a result, in our analysis, we may consider 
5 to be always positive. 

Next, we will try to determine all the singular solutions of the theory with the singu- 
larity occurring at finite time hoping that they do not cover the whole phase space leaving 
some room for the non-singular ones. The whole treatment will be analytical assuming a 
polynomial dependence of the coupling function ^(0) on the scalar field, ^(0) = 0" with n 
being a positive integer greater than unit. Since a singular solution is characterized by the 
vanishing of the scale factor at some finite time, a{t) = e'^*-*-' 0, we will always demand 
that near the singularity y = e"^*^ — > 00. In the same region, the quantity z = u will be 
set to approach a constant value, zero or infinity while the scalar field will be left free to 
adopt any possible behaviour. 

We are going to concentrate our attention on the study of the following cases : 



(I) z = finite 7^ 0, = any, y 00. 



If, near the singularity, z remains finite adopting a constant value, z{t) = c, then, from 



the differential equation for y ( |1.11| ), we obtain 

-2z ^ y = ^ e-^I ^W'^* R{t) ~ e' 



y 
y 



ct 



(3.4) 



which goes to zero only when ct — > —00. This means that the singularity is approached only 
at infinite time and for this reason it must be excluded. 
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(II) z ^ 0, (j) = any, y ^ oo. 

In this case, the first derivative of the scalar field with respect to time takes the form 



3 / 3 

X = — - 6^'zky + Y (2 ^i'^^^uY + 6% . (3.5) 

We have to consider the following cases: 

(A) (| 5i'zkyf > 6% ^ I {6^' zf ky > C»(l). Then, x can be written as 

3 . 35 I >./ 7 I ^sk 

X ^ —-ot,zky + —zyo\t,k\ + 



if s = +1 ^ x+ = , 
= I (3.6) 
if s = — 1 ^ x_ = —36^'zky . 

The above values of x have been taken for the case C,'k > 0. If we change the sign of ^' 
according to ^' — > — we obtain x± — > — Xzp. On the other hand, if we change the sign of k 
in the same way, we are led to x± —>■ ^x^. Note that, apart from the interchange x+ 
the only thing that changes is the absolute sign of x which is not going to be used in the 
following analysis. For this reason, we may consider only the cases ^' > and A; > 0. 

We are going to study each expression of x separately: 
(a) X = X-^- = Then, eq. ( |1.1(]| ) reduces to : 



d^6^'z - I 3{6^'zykyi ■ ^ ' 

8 6^" 

If 3 i^^^'z^ky ^ ^^^^■i th^^' rearranging the differential equations for z and y, we obtain : 

2yz = zy =^ z"^ ^ y which is inconsistent with our assumption for the behaviour of z near 

8 6^" 

the singularity. Next, we assume that „ / r^, \d, — = 6 ^ (9(1). In the same way, we obtain 

3 (^4 z) ky 

~ which is consistent with our assumptions only for (1 — 6) < 0. However, the 

differential equation for y gives the result y^^^^^ ~ (t + c), which leads to the conclusion that 
the singularity is approached only at infinite time and for this reason it must be excluded. 
8 b^" 

If, finally, 3 [^^If ky > ^(1)' are led to 

. 8^r„ ^ ^(^^)_8 m; 



-2yz 
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di^__f y' 
yW 



=> y" m'y - y" - y' y + ^ y = o. (3.8) 

If we assume that C,{(t>) = 0", the only solution of the above differential equation, compatible 
with this assumption, is the following 



provided that 



bi= ^l^"".^ , and bl + 2b2h-h,=0 . (3.10) 



But, then, we obtain : 



{sez? y = m' (-^) y = -y ^e' = -hi = o{i) (3.11) 

which is inconsistent with our assumption that ^ 0{1). 

(b) X = x_ = —36^'zky. Then, eq. ( |1.10| ) takes the form 

i = {-36^'zky) ~ -z^ (1-6 6Cky) . (3.12) 
d(f) 

If we assume that 6 5$,"ky <^ C>{1) or 6 6^"ky = 6 ~ C^(l)) ^6 obtain exactly the same 
result as in case (a). The first assumption leads to infinite z while the second one leads to 
a singularity which is approached only at infinite time. The third option, Q5^"ky ^ 0{1), 
leads to 

dz 

z=— {-36^'zky) ^z^Q dC'ky ^ z- (5^')"' (3.13) 

Z(p 

which goes to zero only if oo. By using the above result and for C,{(f)) = 0" with n > 2, 
the differential equation for y gives : y ~ {SC,")~^ or equivalently that 6^"y ~ (9(1) which is 
inconsistent with our assumption that 6C,"y ^ (9(1) • For the special case of = 2, the same 
equation leads to ?/ ~ Incp and consequently to {S^'z)^y <ti 0{1) which is again inconsistent 
with the assumption {S^'z)^y ^ 

(B) {^6^'zkyY < 6ky ^{6^'z)^ky < C(l). Then, x ~ ±V^, which means that 
the following analysis is valid only for k = +1. The differential equation for z takes the form : 

, = ^,^_,2_ i8-65ey)y .3^. 

del.''- ' 4+|(5av ■ ^ ' 

Now, we have to consider the following cases concerning the quantity 6C,'y that appears 
in the denominator of the above equation : 
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(a) = a = 0{1). This means that — > and consequently that 6^"y ^ 0{1). Then, 
the differential equation for z takes the form 

. ^ 6^ ^dz^±y^ d(5e) y , (3.15) 

where A = 4+ |a^. Since 6^'y = const., we easily obtain d{6C,') y = —6^'dy. Then, we find 
that z ~ which once again leads to a behaviour of z radically different from the assumed. 

(b) 6^'y » 0{1). In this case, and according to our assumption that ^(0) = 0", we can 
only have : S^"y » 0{1). Then, we obtain 



dz 4(5^ 

di±V6^) ^dx 
d^ ^ 



^"^+(^ = 0- (3-16) 



The only solution of the above differential equation, compatible with our assumption for 
i{4>)-, is the following 

bi , ,2 16 (n- 1)2 

X = , where 6? = — '- . 3.17) 

./W ^ n{2-n) ^ ' 

However, the above result leads to 5^"y ~ C'(l) which is inconsistent with our assumption 
that 5S,"y ^ C^(l)- Moreover, the above solution for x is real only if n < 2 which is in 
disagreement with our assumption for ^(0). 

(c) 5^'y -C 0{1). In this case, the differential equation for z takes the simple form z ~ 
— (2 — \ 5^"y)y. For the assumptions 5$^"y <^ and S$,'y = 6 ~ C^(l)i we obtain the 

result y which is different from our assumption that z goes to zero near the singularity. 
The other option, 5^"y 3> C(l), is a little more complicated as it leads to 



d{±yQy) = dx 
d^ 'dd 



x" + ^^0. (3.18) 



In the same way, the only solution of the above differential equation, for ^(0) = 0", is the 
following 



This means that the above solution is real only for n < 1 which is inconsistent with our 
assumption for ^(0). 
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C) {^5S,'zkyf ~ 6% ^ '^{Si'zfky = \ ^ 0{l). Then, x ~ '^Si'zky, where A 
— 1 ± A/r+~a, and eq. ( p. .101) takes the form : 



-z 



o \2 

'[{l + 2a)-—5i"ky} . (3.20) 



Since (5^'z)^?/ ~ 0{1), we will always have 5^"y ^ C'(l) independently of the behaviour of 
the scalar field near the singularity. By using the above in the differential equation for 
we obtain that z ~ {^i'Y which goes to zero for (A > 0, — >• 0) or (A < 0, — > oo). 
The solution of the differential equation for y is y ^ (^^")~^ + c which leads to y ~ C>{1) 
if ^ oo or to yS^" ~ C>{1) if ^ 0. The first of these results is inconsistent with 
the assumed behaviour of y near the singularity while the second one disagrees with our 
assumption that y6C," ^ C'(l)- 

From the study of the first two cases, (I) and (II), we conclude that in a singular 
cosmological solution the first derivative of the quantity u cannot remain finite or vanish 
near the singularity. The only option left is the adoption of an infinite value, just like y, 
which we are going to study now. 



(Ill) z oo, = any, y ^ oo. 

We are going to study the following cases : 

(A) (| 5ezY {z^ + kyf > 6 (z2 + ky) (SCz)^ {z^ + ky) > C(l). 

As in case (II), the change of sign of ^' or z'^ + ky leads to the interchange x+ ^ x_ as well 
as to the change of the absolute sign of x. However, the sign of x does not appear anywhere 
in our analysis while the interchange between the two expressions of x does not affect our 
arguments since both of them are being studied. As a result, we are going to assume again 
that both of ^' and z'^ + ky always take on positive values. Then, under the above assumption, 
the quantity x assumes the expressions x+ = 2/{6^'z) and X- = —36^z {z^ + ky) for s = ±1, 
respectively. 



(a) x = x+ = Then, 



dz 2 2 f. 8 



5i'z 1 3 {di'zf (z2 + ky) 



(3.21) 



First, we assume that , .a /" >) ; — r -C Oil). Rearranging the system of differential 

{o^ z)^ {z^ + ky) ^ ' 

equations for z and y, we easily obtain that z^ y which means that the two quantities 
have exactly the same behaviour near the singularity and, if necessary, their sum may be 
written as z^ + ky = az"^, where a an arbitrary constant. The differential equation for z gives 
the result z^ = {6^ + c)~^ which goes to infinity when (5^ + c) — > 0. However, we have to 
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check that this hmit takes place at finite time. From the expression of we obtain 



-^ = j^^JS^ + c^ {5i + c)^{t + df (3.22) 

which indeed goes to zero at finite time. It is easy to check further that our basic assumptions 
{z^ + ky) > 0(1) and (^^^/^)/(^^ ^ j^^) < ^(1) are satisfied provided that 5^' ^ 0. 
As a result, we conclude that the solution 

is an acceptable singular cosmological solution with the singularity being approached at finite 
time. 

Next, we consider the case / ^y., / '> ; — \ = b ~ Oil). Since both z and y are 

[di zY [z^ + ky) ^ ^ 

infinite near the singularity and in order to satisfy the above constraint, we must have 
S^' 0. The corresponding solution for z is z~'^ = (1 — 6) {S^ + c). If 2;^ <^ y, we obtain 
{6^'zf {z^ + ky) ~ {6^'z)^y^ 0{1) and that z^ - y^^-^\ Then : 



{5^'z)^y 8 



{6a~'^z'y + c' . (3.24) 



The above result leads to {6^'z)'^y ~ 0{1) which is inconsistent with our assumption. Fol- 
lowing an exactly similar analysis, we may show that when z^ >> y or ~ y, the above 
integration yields a constraint which is again in obvious disagreement with our assumption 
that {S^'z)^z^ + ky) » C(l). 

Finally, we assume that ^^^/^y ^^-2 _|_ ^ ^(1)- furthermore, z'^ ~ y, we may set 

z'^ + ky = az"^ and the differential equation for z gives z^ ~ ((^^')~^- ^'^^i this leads to 
{5^' z)"^ [z"^ + ky) ~ 0{1) which is inconsistent with our assumption. If z"^ ^ we reach the 
same result as in the case z"^ ~ y. Finally, the case z"^ <^y has been studied in (IIAa) and 
has been shown to lead to the result {5^' z)"^ {z"^ + ky) ~ 0(1) which is inconsistent with our 
assumption. 

(b) X = x_ = —"iS^'z {z"^ + ky). In this case, the equation for z takes the form 

Az^ + iz^ + ky)j ■ ^"^-^^^ 

We start by assuming that ^^2 _|_ ^^2 _|_ ^ 0{1). Then, from the differential equations 

for z and y, we obtain z"^ y which leads to S$," z^ -C C(l). Setting z^ + ky = az^ and 
integrating the equation for z, we find that SC,' z"^ ~ + ... from which we conclude that 
the scalar field must go to infinity near the singularity. Unfortunately, this result leads to 
5C"z'^ ~ 5C'(f>{^C)'^ ^ 0{l) which is inconsistent with 6^' z^ < 0(1). 
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If we assume that \ = 6 — 0(1) and, furthermore, that z'^ ~ y, the 

differential equation for z gives the result : 5^' ... which goes to infinity only when 

2a (1 - h) 

(f) ^ 00. Substituting this in the same equation, we obtain that z ~ ^(^^ + 46 — Sab) 
order to assure the assumed behaviour of z near the singularity, we demand the positivity 
of the exponent of (f). But, then we conclude that ^'(0) < (f) which is inconsistent with our 
assumption that ^{(f)) = 0" with n > 1. For the case z"^ ^ y, wc obtain the same result with 
a = 1. Finally, if z"^ ^ y, the two constraints become 6^'y ^ 0{1) and QS^"y = b. From 
the differential equation for y, we obtain that y6^' ~ (p + ... which again goes to infinity 
only if — > 00. Substituting this result in the same equation, we find that y ~ with 
4 + 6 > 0. But, this leads once again to the inconsistent result ^'(0) < 0. 

Finally, we consider the case where g^^^^f^ » 0(1). If z'^ ~ y, the differential 
equation for z gives z ~ (5^')-V(4+a) ^y^j^ j^g^^g ^^/^2 _ (^^^/)(4-3a)/(4+a)_ jf g^^^ 
y = (a — l)z'^, the differential equation for y gives the result z ~ 02a/(3a-4)_ jf ^ > 4^/3 
and 0^0, the constraint t^^'z^ » C'(l) is satisfied but 2; — > 0. The same holds if a < 4/3 
and (p —>■ 00. In both cases, the behaviour of z is inconsistent with the one that has been 
assumed. If z"^ ^ y, the equation for z gives z ~ {S^')'"^^^ which goes to infinity only when 
S^' — > 0. But then, S^'z"^ ~ i^CY^^ ~^ which does not agree with our assumption. Finally, 
if <^ y, we obtain that z ~ {S^')^^ which goes to infinity only if — > 0. Then, the differ- 
ential equation for y gives y ~ + c or equivalently yS^" ~ 0{1) which is inconsistent 
with our assumption. 

(B) (| S^'zf {z^ + kyf > 6 + ky) {S^'zf {z^ + ky) < 0(1). 



In this case, the expression of x becomes x = ±^^6 {z'^ + ky) while the differential equa- 
tion for z takes the form 




- 6^^^^ {z^ + ky) ] {z^ + ky) 
4z^ + f (S^'z)^ {z^ + ky)^ 



(3.26) 



In order to simplify the analysis, we are going to consider first the possible relation 
between z'^ and y and then study each case separately. 

(a) z^ ~ y. Then, the constraint becomes S^'z"^ -C 0{1) and the differential equation 
for z reduces to 

(l + 2a)--a2 6ez^ 



2 

2; ~ — 2; 



2 



(3.27) 



where we have set z^ + ky = az^. If we assume that S^" -C C(l), the above equation gives 
the result z ~ which goes to infinity only if — » ±00. But, then, S^' 00 and 

there is no way that the constraint 6^'z^ ^ C^(l) can be satisfied. A similar result arises 
when I 6i"z'^ = 6 ~ 0(1). If we assume that 5i"z'^ > 0(1), we find that z"^ ~ {5^' + c)"^ 
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which goes to infinity only when 6^' + c ^ 0. But, then, the quantity S^'z"^ goes to infinity 
as well and this is inconsistent with our initial assumption. 

(b) z'^ ^ y. In this case, the results are similar to the ones of the previous case and 
they easily arise from them if we set a = 1 . 

(c) z"^ <^ y. In this case, the constraint becomes {S^'zYy <S 0{1) and the differential 
equation for z takes the form 

2 (8-65^"ky)ky 
zc^-z^- \ , s,rJ\o ■ 3.28 

The above equation has been arisen also and studied in case (IIB). There, it was shown 
that, for the cases 6^'y 3> or -C there was no real solution of the corresponding 

differential equation for x compatible with the assumption that ^(0) = 0" with n > 2. 
On the other hand, for the case S^'y = a ~ 0{1), we were led to 2; ~ which is again 
inconsistent with our assumption that z"^ <^y. 



(C) (| Si'zf {z^ + kyf ~ 6 + ky) ^ | {Si'zf {z^ + ky) ^ ^ 0{1). 

Since both of z and y approach infinity near the singularity, the only way to fulfill the 
above condition is to have 0^0. The expression of x becomes x — ^ ^S^'z {z"^ + ky), where 
A = — 1 ± Vl + o, while the differential equation for z takes the form 

. _ 2 _ (8a - d^'ax^) {z^ + ky) 

^--^ - Aa-8X + A{z'' + ky)/z^ ' ^ ^ ^ 



We consider the cases : 

(a) z"^ y . Then, we may set z"^ + ky — Pz"^ — Py and the equation for z reduces to 



We assume further that 5^"z'^ < 0{1) or equivalently that 5Cy < 0{1). Then, the 
differential equation for y gives the result S^'y ~ (0~^ + c) which leads to 5^'y ~ 0{1) or, 
since y ~ z^, to S^'z"^ ~ C^(l)- Substituting this in the above equation for we obtain that 
z ~ e~^'^, where i? is a combination of a, A and j3. This means that near the singularity, 
2; — > which is inconsistent with the assumed behaviour of z. If we assume that 6^"z'^ = 6 ~ 
(9(1), we are led to exactly the same result. Finally, we consider the case S^"z'^ » C'(l). The 
equation for z leads to the result z ~ (^(5^')2A/3/(2a-4A+2/3) order to fulfill the condition 

that S^'z"^ ^ C^(l) we demand that 2a — 4A + 2/3 + 4A/9 = 0. By using the dependence of 
A on a, we find that the solution of the above algebraic equation is : (3 = 1 = —a. This 
means that z"^ + ky = Pz"^ = z"^ or equivalently that z"^ ^ y which is inconsistent with our 
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assumption. 

(b) ^ y . Now, the constraint becomes {5^' z'^Y — while the differential equation 
for z is the same as in the previous case where we have set /5 = 1. The first two options, 
namely 5^" z^ -C and 5^"z'^ = 6 ~ 0{1), lead to the inconsistent result z ^ The 
third option, 5^" z^ » C'(l), leads, as above, to an algebraic equation with the sole solution 
a = — 1. However, from our constraint, it follows that a is always a positive constant which 
reveals the inconsistency of this solution as well. 

(c) z'^ <^y . Then, the differential equation for z takes the form 

zc^-z'l^{l + 2a)-^5ey^ . (3.31) 

We start by assuming that 6^"y -C 0{1). In this case, the differential equation for y gives 
the result S^'y ~ + c ~ C(l) since, as we concluded above, (p —>■ near the singularity. 
Then, we obtain that (S^'z^y ~ (z'^/v) <S ^{^) which is clearly in disagreement with our 
initial assumption. The same result arises for the case S^"y = b ~ Finally, we assume 

that SC,"y ^ Then, as in case (IIC), we obtain that y ~ {S^")^^ + c or equivalently 

that SC,"y — 0{1) which is again inconsistent with our assumption that S^"y ^ 

So, the only singular solution found, with the singularity occurring at finite time, has the 
form 

where 0^ and tg stand for the value of the scalar field and time, respectively, at the singularity. 
Since y = e~^'^, the scale factor of the universe in this case behaves near the singularity as 

a{t) (t-ts) — > when t ts (3.33) 

which corresponds to a linearly expanding universe with a true cosmological singularity at 
t — ts- During our analysis, no restriction on the value of the parameter k has been arisen 
and, for this reason, the above singularity can be considered either as an initial or as a 
final one. Note, that the rate of expansion is much larger that the corresponding ones 
during the "radiation" or the "matter" epoch of the Standard Cosmological model. The 
same result was found in the first section in the case of an empty universe. However, when 
the presence of the energy momentum tensor of the scalar field was taken into account the 
linear dependence of the scale factor on time t changed into the milder dependence ~ t^^^. 
Now, after the addition of the Gauss-Bonnet term in the theory, we find again a scale factor 
linearly dependent on time. This could be interpreted in the following way : although the 
addition of the Gauss-Bonnet term has made the search for singular cosmological solutions 
much more complicated, the net result of its presence in the theory is to cancel exactly the 
contribution of the energy momentum tensor of the scalar field leaving behind an "empty" , 
linearly expanding universe. The only difference between these two cases is that while the 
truly empty universe gives rise only to an open universe, as one would expect, the virtually 
"empty" Gauss-Bonnet universe can be interpreted as a closed or open universe. 
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The question which arises next is whether the presence of the Gauss-Bonnet term affects 
the singular cosmological solutions in such a way that they are restricted in only a region 
of the phase space of the theory leaving some room for the non-singular ones. The singular 
solutions ( |3.32|) were determined under the assumptions x = and 6^' 7^ 0. Then, x 

tends to zero only asymptotically, when the singularity is approached, keeping a definite 
sign otherwise. If we assume that the singular solutions are characterized by x > 0, then 
we must always have ^' > since 6 is always positive. If C,{4>) = 0" with n even, the first 
derivative and hence x, is positive only if is positive as well. This leads to the restriction 
of the singular solutions in the positive (p half-plane in a z-(f) graph leaving the other half for 
the non-singular ones. The above argument leads to the selection of even values of n since 
for odd values the first derivative of the coupling function retains always a positive sign, for 
all values of 0, and this leads to the non-confinement of the singular solutions in a part of 
the phase space of the theory. 

Finally, it is worth noting that for the cases k = and k = —1 the singular solutions are 
restricted not only in the positive half-plane for x > and 6 > but, moreover, in the 



positive z quarter-plane. This can be justified if we examine eq. (|1.9| ) written in the form 

{z^ + ky) (6 - ?>5i'zx) = x^ . (3.34) 

If we assume that k = +1, then, the value 2; = is indeed acceptable by the above equation 
which means that, in this case, the singular solutions can traverse the axis of at the point 
where = 6y and extend at both the upper and lower z quarter-planes. On the other 
hand, if = — 1, the value z = must be excluded in order to ensure the reality of x. This 
means that, now, the singular solutions remain confined in one of the two z quarter-planes 
according to the sign of z that has been initially chosen. Finally, if = 0, setting 2; = in 
the above equation leads to the result that the singular solutions traverse the 0-axis at the 
point where x = 0. However, this point is only reached near the singularity where 2; ^ 00. 
As a result, the value z = is not an acceptable one for the option k = and the singular 
solutions are again restricted in one of the two z quarter-planes. 



4 Numerical Analysis 

In this section, we are going to demonstrate numerically the existence of non-singular 
solutions in the presence of the Gauss-Bonnet term in the action functional of the theory. 
The analytical study of the previous section has led to the result that the singular solutions, 
characterized by the occurring of a cosmological singularity at finite time, exist in the context 
of the theory but they do not cover the whole phase space. As we will see, it is the non- 
singular solutions of the theory that are summoned to fill the rest of the space. 

As we have already mentioned, the coupling function of the scalar field to the quadratic 
Gauss-Bonnet term is of the form ^ (0) = 0" with n even. As an illuminating example, we are 
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going to assume the simpler case, that is n = 2, in our numerical analysis. The system of the 
differential equations ( |1.10|) -( p!.ll| ) together with the algebraic equation (|1.12| ) is numerically 



integrated yielding a solution for z = uj and a = . Since the case A; = has already been 
studied in Ref. |]^, we are going to concentrate our attention on the other two cases, namely 
k = ±1. 

The solution for the case = +1 is depicted in Figures 1 and 2 for a family of singular and 
non-singular solutions. As we have analytically proven in the previous section, the singular 
solutions, in a z-(f) graph, are confined in the positive half-plane, with the ^-axis playing 
the role of a barricade for them, since we have chosen 6 = 0.5 > 0. Changing gradually 
the boundary conditions of the numerical integration, the singular solutions cease to exist 
as we approach the z-axis while the non-singular solutions start to develop and cover the 
negative half-plane. According to our analysis of section 3, the singular solutions are 
characterized by the simultaneous divergence of the quantity z and of the quantity y or, 
equivalently, the vanishing of the scale factor a. As we also noted, both of the singular 
and non-singular solutions can traverse the axis and extend at the upper as well as at 
the lower z quarter-plane. The value z = 0, for k = +1, is an acceptable value for both 
families since they obey the same equation ( p.34| ). The dependence of the scale factor of 



the universe is displayed in Figure 2 where the avoidance of the initial as well as of the final 
singularity of the closed universe is obvious. Both of Figures 1 and 2 have been drawn for 
the choice s = +1. According to our argument of section 2, the solutions, singular or not, 
of the equations ([LToD - dTTlD are invariant under the simultaneous change of z and s. This 



means that if we choose s = — 1, the Figures 1 and 2 will remain unchanged apart from the 
sign of z in Figure 1. 

The asymptotic form of the non-singular solutions, for the case k = +1, for early and 
late times can be found by making use of the ansatz 

= 00 + 01^'' , CO = CUo + CUit" . (4.1) 



Substituting the above expressions in the system of equations ( |1.4| )-(]L^), we find that the 
only acceptable non-singular solution has the form 

= 00 + ^ , uj=-ln{ — )+^, (4.2) 



VS ' 2 \2 

where 0o and loi are arbitrary constants. According to the above solution, the asymptotic 
regions — *• ±oo correspond exactly to t — >• ±oo modulo a constant coefficient. Then, as it 
is obvious from Fig. 2 as well, the non-singular solutions are characterized in both limits by 
the same constant u parametrized in terms of the parameter S, or equivalently by the same 
constant scale factor a, which means that our solutions interpolate between the same static 
Einstein universe without undergoing a cosmological singularity, either initial of final. 

A family of non-singular solutions, for the choice k = —1, is depicted in Figures 3 and 4. 
As we note in Figure 3, the singularity-free solutions cover the whole positive z half-plane 
without being able to cross the 0-axis since this would lead to an imaginary value of x 
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according to eq. ( |3.34| ). The singular solutions, although not shown in Figure 3, cover the 
lower-right-hand-side quarter-plane of the graph, for the choice 6 = 0.5 > 0, and they are 
completely separated from the non-singular ones. In Figure 4, the dependence of the scale 
factor a on the scalar field is displayed with the absence of the initial singularity of the open 
universe being obvious. The Figures 3 and 4 have been drawn for the choice s = — 1. As in 
the previous case, the alternative choice s = +1 leads to the same graphs apart from the sign 
of z due to the corresponding invariance of the equations ( p..lO| )-( p!.ll| ). According to Fig. 4, 
the scale factor adopts a non- vanishing, constant value in the limit (p oo while it increases 
rapidly as — oo. Given that the dependence of the scalar field (p on time t is a smooth 
function leading to the asymptotic behaviour — > ±cxo when t — > =Foo, we may conjecture 
that the non-singular open universe interpolates between a static Einstein universe at early 
times and an expanding universe at late times. 

In order to complete the picture, we display the graph of the gravitational scalar curvature 
R = gf^'^ R^y = 6 ci) + 12 cj^ + Qke""^^ versus the scalar field 0, for the cases of the closed and 
open universe, in Figs. 5 and 6, respectively. The absence of divergences, which implies 
the absence of cosmological singularities, is obvious in both figures. For the case of the 
closed universe, k = +1, the scalar curvature interpolates between two constant values, as 
(j) — > ±cxD, while in the case of the open universe, k = —1, the scalar curvature adopts a 
constant value, in the limit (p (yo, reaches a maximum value near = and vanishes 
quickly as we approach the asymptotic region —oo. 



5 Conclusions 

In this article, we have studied a general field theory that describes the coupling of a scalar 
field to higher-curvature gravity through the quadratic Gauss-Bonnet term. As we men- 
tioned in the Introduction, the presence of this term to the action functional of the Super- 
string Effective Theory has led to the existence of new black hole solutions and non-singular 
cosmological solutions. For reasons that will be clarified shortly, we have chosen to study 
a slightly different quadratic gravitational theory where the coupling function between the 
scalar field and the Gauss-Bonnet term has the polynomial form ^(0) = 0", with n being a 
positive integer. 

The classical, scalar and gravitational, equations of motion were solved initially in the 
absence of the Gauss-Bonnet term. Assuming a constant or a time-evolving scalar field, 
we determined the cosmological solutions that correspond to an empty universe or to a 
universe that is dominated by the energy-momentum tensor of a scalar field, respectively. In 
both cases, the cosmological solutions were found to be characterized by, at least, one true 
cosmological singularity which could not be avoided if a realistic solution was demanded. 

In the presence of the Gauss-Bonnet term, the energy conditions were shown to be 
violated allowing the existence of non-singular cosmological solutions. For k = 0, +1 this 
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violation holds only for 6 > while for k = —1 the conditions are violated for both signs of 6. 
However, only for 6 > non-singular solutions arise in the context of the theory. Next, the 
equations of motion were solved analytically near the region of the cosmological singularity. 
A family of singular cosmological solutions was determined and was shown to correspond to 
a linearly expanding universe, a(t) ~ (t — ts), with a true cosmological singularity at finite 
time ts- This singularity can be considered as an initial or a final one and the corresponding 
universe can be interpreted as an open or closed. One of the basic conclusions drawn from 
this analysis was that these solutions are confined in a part of the phase space of the theory 
leaving the rest of the space for the non-singular ones. This remarkable feature holds only 
when the coupling function has the dual, polynomial form ^(0) = 0", where n is an even, 
integer number. 

The specific choice ^ (0) = 0^ was studied numerically and the solution for the parameter 
z and the scale factor a was determined. For both values of the parameter k, the singular 
cosmological solutions were confined in a part of the phase space of the theory, as predicted, 
and completely separated from the non-singular ones. For the case k = +1, was even 
possible to found the asymptotic, analytic form of the non-singular solutions in the limit 
t —>■ ±oo. According to this form and our numerical results, the universe, for the case 
k = +1, interpolates between the same asymptotic, static Einstein universe with constant 
scale factor passing through an intermediate phase of expansion-contraction. By making 
use of the asymptotic form of the non-singular solutions ([4.2|) for early and late times, we 

may easily find that the scale factor asymptotically adopts the value : a{t) = e"^ ~ W^. 
This means that the coupling parameter 6 plays the role of the asymptotic value of the scale 
factor of the universe which explains the appearance of non-singular cosmological solutions 
only for positive values of 6. According to the above, in the limit 5 — >■ 0, that is when we 
eliminate the Gauss-Bonnet term from the action functional of the theory, the asymptotic 
value of the scale factor goes to zero and we recover the singular cosmological solutions of 
section 2. For the case k = —1, we may conjecture that the universe interpolates between 
a static Einstein universe at early times and an expanding universe at late times. For both 
values of k, the absence of cosmological singularities, a feature that owes its existence to the 
presence of the higher-derivative Gauss-Bonnet term, is obvious. 

As is well known, similar results arise in the context of the Superstring Effective Theory 
where the coupling function between the scalar field and the Gauss-Bonnet term has the 
form 

^{(f)) = ln[2e'^T]^{ie'^)] , (5.3) 

where t] is the Dedekind function. However, the above form and the polynomial form that 
we have chosen during this article share a number of important characteristics. More analy- 
tically : they are both invariant under the change —0, they both have a global minimum 
at = and they both take on an infinite value as — *■ ±oo. As a result, we may 
state the following "theorem" : Any theory of the form that describes the coupling 

of a scalar field with the Gauss-Bonnet term through a coupling function which (i) is dual, 
(a) has a global minimum and (Hi) asymptotically tends to infinity may lead to singularity- 
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free cosmological solutions. Superstring Effective Theory is a characteristic example of such a 
theory but it is not the only one. Another choice, as we have demonstrated both analytically 
and numerically, is the theory ( |1 . 11 ) with ^{(p) = 0^" and a number of alternative choices 
may follow as long as the corresponding coupling functions satisfy the three aforementioned 
criteria. 

An important subject which has not been addressed in this article is the stability of our 
singularity-free solutions. Here, we have focused our attention to the existence itself of these 
solutions. The question of their stability under linear or even non-linear perturbations is 
an independent subject which still remains open. If our singularity-free solutions turn out 
to be stable, that will constitute a great achievement towards the resolution of the initial 
singularity problem. If, on the other hand, these solutions turn out to be unstable, their 



collapse in the early universe may have led to the creation of primordial black holes [ITU] . We 
hope to come back to these open questions in a future publication. 
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Figure 1: The dependence of the quantity z on the scalar field for a family of singular and 
non-singular solutions for the case A; = -|-l, s = and (5 = 0.5. 
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Figure 2: The dependence of the scale factor a of the universe on the scalar field for a 
family of singular and non-singular solutions for the case A; = -|-l,s = -|-l and 5 = 0.5. 
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Figure 3: The dependence of the quantity z on the scalar field for a family of non-singular 
solutions for the case k — —1, s — —\ and b = 0.5. 
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Figure 5: The scalar curvature R versus the scalar field for a family of non-singular 
solutions for the case k = +1, s = +1 and 6 = 0.5. 
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Figure 6: The scalar curvature R versus the scalar field (j) for a family of non-singular 
solutions for the case k = —1, s = — 1 and S = 0.5. 
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